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Abstract. - Monte Carlo (MC) and series expansion (SE) data for the energy, specific heat, magne- 
tization and susceptibility of the two-dimensional 4-state Potts model in the vicinity of the critical 
point are analysed. The role of logarithmic corrections is discussed and an approach is proposed 
in order to account numerically for these corrections in the determination of critical amplitudes. 
Accurate estimates of universal amplitude ratios A+/A-, r+/r_, r-r/r_ and are given, which 
arouse new questions with respect to previous works. 



Introduction. The concept of universality is of fundamental importance in the theory of 
phase transitions. Critical exponents and critical amplitudes describe the leading singulari- 
ties of physical quantities in the vicinity of the critical point, 

M_(r) « B(-rf, X± (r) » r ± |r|^, C±(r) « ^\r\~ a , (1) 

a 

(r = (T — T c )/T is the reduced temperature and the labels ± refer to the high-temperature 
(HT) and low-temperature (LT) sides of the critical temperature T c ) and universal combi- 
nations of critical amplitudes [1], as well as critical exponents characterize the universality 
class of the model. For the Potts models with q > 2, in addition to the above mentioned 
quantities, a transverse susceptibility is defined in the LT phase. 

Analytical results for the critical amplitudes for the g-state Potts models with q = 1, 2, 
3, and 4 were obtained by Delfino and Cardy [2], using the two-dimensional scattering field 
theory of Chim and Zamolodchikov [3] . In the case of the 4-state Potts model, the approach 
of Rcf. [2] leads to the universal susceptibility amplitude ratios T + /F_ = 4.013 and Ft /T_ = 
0.129. Monte Carlo (MC) simulations also reported in [4] did not confirm conclusively these 
predictions. Another MC study due to Caselle et al. [5] leads to r+/T_ = 3.14(70), which is 
below the theoretical prediction of Delfino and Cardy. More recently Enting and Guttmann 
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analysed new (longer) series expansions for q = 3 and q = 4 obtained by the finite lattice 
method [6]. Their estimates T + /r_ = 3.5(4) and r T /T_ = O.ff (4) for q = 4 arc in slightly 
better agreement with the results of [2] and [4]. An analysis by differential approximants 
however is successful only for q = 3 where the corrections to scaling are represented by pure 
powers, but meets with some difficulty in the q — 4 case, in which logarithmic corrections 
are expected. Therefore they had to resort to a slowly convergent direct analysis of the 
asymptotic behaviour of the expansion coefficients. 

In this letter we present accurate Monte Carlo data supplemented by a re-analysis of the 
extended series derived in [6] . We are essentially concerned with the universal combinations 

i± E± El (2) 
a_' r_' r_' [ ' 

To the various critical amplitudes of interest A± , . . .we have associated appropriately defined 
" effective amplitudes" , namely temperature dependent quantities A± (r) , . . . which take as 
limiting values, when |t| — > the critical amplitudes A±, .... By analogy we have also con- 
sidered " effective ratios" of critical amplitudes which tend to universal ratios as r — > and 
exhibit smoother behaviours in the vicinity of the critical temperature than the quantities 
themselves. Considering effective ratios would even eliminate logarithmic corrections from 
the fit in the case of 4-state Potts model in absence of regular contributions, which unfor- 
tunately do exist! We also use the self-duality relation to check explicitly the cancellation 
of the dominant corrections to scaling in the case of the energy density evaluated at dual 
temperatures. 

Model and observables. The Hamiltonian of the Potts model reads as H = — ^sisj , 
where takes integer values between and q — I, and the sum is restricted to the nearest- 
neighbor sites (ij) on the lattice. The partition function is defined by Z = J^conf e ~^ H ■ Cm 
the square lattice, in zero field, the model is self-dual. The duality relation (eP — l) (eP — f) = 
q determines the critical value of the inverse temperature (3 C = ln(l + ^/q). Dual reduced 
temperatures r and r* can be defined by (3 = /3 C (1 — r) and (3* = /3 C (1 + r*). 

In our simulations we use the Wolff algorithm [7] for square lattices of linear sizes L = 20, 
40, 60, 80, 100 and 200 with periodic boundary conditions. Starting from an ordered state, 
we let the system equilibrate in I0 5 steps measured by the number of flipped clusters. The 
averages are computed over 10 6 — 10 7 steps. The data are measured in a range of reduced 
temperatures called the "critical window" and defined as follows: the lower limit is reached 
when 1 7" | ^ reaches the size L of the system, and the upper limit of the critical window is 
fixed for convenience when the corrections to scaling in the Wegner asymptotic expansion [8] 
do not exceed a few percent, say 2 — 3%, of the leading critical behaviour Eq. (I) (forgetting 
about the logs). This definition avoids finite size effects which would otherwise make our 
analysis more complex. 

The order parameter of a microstate M(t) is evaluated at the time t of the simulation as 
M = qNn ^l~ X , where N m is the number of sites i with si = m and m £ [0, q — 1] is the 
spin value of the majority state. N = L 2 is the total number of spins. The thermal average 
is denoted M = (M). Thus, the longitudinal susceptibility in the LT phase is measured by 
the fluctuation of the majority spin orientation ksTx- = (N^) — (N m ) 2 and the transverse 
susceptibility is defined in the LT phase in terms of the fluctuations of the minority of the 
spins k B TxT = j^jj E^m((^ 2 ) ~ (^) 2 )> whilc in thc HT P hasc X+ is given by the 

fluctuations in all q states, ksTx+ = ^Ej=o((^) ~ (^/j) 2 )i where is the number of 
sites with the spin in the state [i. This definition of the susceptibility is, in both phases, 
completely consistent with the available series expansion data [9]. The internal energy 
density of a microstate is calculated as E = —jj'^(ij) ^sis,, and its ensemble average is 
denoted as E = (E). The specific heat measures the energy fluctuations, (fcsT) 2 C = — ^§ = 
«E 2 )-(E) 2 ). 
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Our MC study of the critical amplitudes will be supplemented by a rcanalysis of the HT 
and LT expansions recently calculated through remarkably high orders by Enting. Guttmann 
and coworkers [6, 10]. In terms of these series, we can compute the effective critical ampli- 
tudes for the susceptibilities, the specific heat and the magnetization and extrapolate them 
by the current resummation techniques, namely simple Pade approximants (PA) and dif- 
ferential approximants (DA) properly biased with the exactly known critical temperatures 
and critical exponents. The LT expansion, expressed in terms of the variable z = cxp(— /?), 
extends through z 43 in the case of E , z 59 in the case of X- (z 47 f° r Xt) and z 43 for M. The 
HT expansions, computed in terms of the variable v = (1 — z)/(l + (q — l)z), extend to v 43 
in the case of E, and v 2i for the x+- As a general remark on our series analysis, we may 
point out that the accuracy of the amplitude estimates is questionable, since the mentioned 
resummation methods cannot reproduce the expected logarithmic corrections to scaling and 
therefore the extrapolations to the critical point are uncertain. In this case we have also 
tested a somewhat unconventional use of DA's: in computing the effective amplitudes, we 
only retain DA estimates outside some small vicinity of the critical point, where they appear 
to be stable and reliable. Finally we perform the extrapolations by fitting these data to an 
asymptotic form which includes logarithmic corrections. 

Logarithmic corrections. In the usual parametrization cos(7rj//2) = \^fq in terms of 
which the scaling dimensions are known, we have y = at q = 4 and the second thermal 
exponent [11,12] = — 4y/3(l — y) vanishes. Accordingly, the leading power-behaviour of 
the magnetization (and of other physical quantities) is modified [13] by a logarithmic factor 

M_(-|t|) = J B|r| 1 /i2 ( _ ln | T | ) -i/8^ corr( _ ln | T | ); (3) 

and a correction function T corr {— In |r|) contains terms with integer powers of (— ln|r|), 
ln(— In |t|)/ (— In |r|),. . . Non-integer power corrections may also occur due to the higher (ir- 
relevant) thermal exponents [11,12,14,15] y < p n or to other irrelevant fields, but let us first 
discuss the form of the logarithmic terms. Extending the pioneering works of Cardy, Nauen- 
berg and Scalapino (CNS) [13,16], Salas and Sokal (SS) [17] obtained a slowly convergent 
expansion of T corr (— In |r|) in logs, e.g. for the magnetization: 



il/_(-|r|) = J B|r| 1 /i2 ( _ ln | T | ) -i/8 



1 3 M-lnM, 



16 -lnM Vln 



(4) 



We provide below a re-examination of this and similar quantities. The non-linear RG 
equation for the relevant thermal and magnetic fields <f> and h, with corresponding RG 
eigenvalues y$ and yh, and the marginal dilution field ?/>, are given by 

{y<t> + y^)<P, (5) 



din b 

j^t = (Vh + Vh^)h, (6) 
dtp 



din 6 



= <?W0 (7) 



where b is the length rescaling factor and / = In b. The function g(ip) may be Taylor 
expanded, g(ip) = y^ip 2 ^ + ^-ip + ...). Accounting for marginality of the dilution field, 

there is no linear term at q = 4. The first term has been considered by Nauenberg and 
Scalapino [16], and later by Cardy, Nauenberg and Scalapino [13]. The second term was 
introduced by Salas and Sokal [17]. For convenience, we slightly change the notations of SS, 
denoting by yij the coupling coefficients between the scaling fields i and j. These parameters 
take the values y^ = 3/(47r), y^, = l/(167r), y^2 = 1/ir and y^,3 = — l/(2-7r 2 ) [17], while 
the relevant scaling dimensions are y<f, = 3/2 and yu = 15/8. 

The fixed point is at <j> = h = 0. Starting from initial conditions </>o, ft-Oj the relevant 
fields grow exponentially with I. The field <p is analytically related to the temperature, so 
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the temperature behaviour follows from the renormalization flow from 6 ~ |t| up to some 
6 = 0(1) outside the critical region. Notice also that the marginal field ip remains of order 
O(ipo) and ipo is negative, |^o| = 0(1). In zero magnetic field, under a change of length 
scale, the singular part of the free energy density transforms according to 

f(^o,M = e- Dl f(iP 7 <f>) 7 (8) 

where D = 2 is the space dimension. Solving Eq. (5) and (7) leads to 

/ = -— lnx+^Mnz, (9) 

where z = ^ v ^+ v ^ and x = d> 16 (for brevity we will denote v = 1/ys = f , u = = 

|) and we deduce the following behaviour for the free energy density in zero magnetic field 
in terms of the thermal and dilution fields, 

f(6 ,i> ) = x Dv z- D »f(6,i>). (10) 

The other thermodynamic properties follow from derivatives with respect to the scaling 
fields, e.g. E(6 ,ip ) = j^f(ip 0l 6 ) = x Dv ^ 1 z~ D > 1 E((j) 1 T/j). What appears extremely 
useful is that the dependence on the quantity z cancels (due to the scaling relations among 
the critical exponents) in appropriate effective ratios. This quantity z is precisely the only 
one which contains the log terms in the 4-state Potts model, and thus we may infer that not 
only the leading log terms, but all the log terms hidden in the dependence on the marginal 
dilution field disappear in the conveniently defined effective ratios. Now we proceed by 
iterations of Eq. (9) and eventually we get for the full correction to scaling variable the 
heavy expression z = const x (— In |r|) £(— In |r|) T(— In |r|), where £(— In |r|) is a universal 
function 



„, , , .. / 31n(-ln|r|)\ / 3 ln(- In |r|)\ 1 / 3 1 \ , , 

while T(— ln|r|) is a function of the variable (— ln|r|) only, where non universality enters 
through the constant ipo- Remember here that x ~ \t\. 

In a given range of values of the t, the function J-(— In |r|) should be fixed and the only 
freedom is to include background terms and possibly additive corrections to scaling coming 
from irrelevant scaling fields. Among the additive correction terms, we may have those of 
the thermal sector = —vy$ n , where the RG eigenvalues are y^ n = D — \n 2 , n 
1.2..'i... [12]. The first dimension = y^ = 3/2 is the temperature RG eigenvalue. The 
next one is y^ = and this leads to the appearance of the logarithmic corrections, such that 
the first Wegner irrelevant correction to scaling in the thermal sector is A^ 3 = —vy^ =5/3. 
In the magnetic sector, the RG eigenvalues are given by yh n = D — \(2n — l) 2 . The 
first dimension y^ = y^ = 15/8 is the magnetic field RG eigenvalue. The second one 
is still relevant, y/ l2 = 7/8, and it could lead, if admissible by symmetry, to corrections 
generically governed by the difference of relevant eigenvalues (yt ll ~yh 2 )/y<t> = 2/3. The next 
contribution comes from yh 3 = — 9/8 and leads to a Wegner correction-to-scaling exponent 
Af l3 = —isyh 3 =3/4. Eventually, spatial inhomogeneitics of primary fields (higher order 
derivatives) bring the extra possibility of integer correction exponents y n = —n in the 
conformal tower of the identity. The first one of these irrelevant terms corresponds to a 
Wegner exponent Ai = — v(— 1) = 2/3 and it is always present. We may thus possibly 
include the following corrections: |t| 2 / 3 , |t| 3 / 4 , |t| 4 / 3 , |t| 5 / 3 , the first and third ones 
being always present, while the other corrections depend on the symmetry properties of the 
obscrvables. 
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In the Baxter- Wu model, which belongs to the 4— state Potts model universality class 1 , 
the magnetization obeys the asymptotic form [19,20] M_{— |r|) = i?|r| 1 / 12 (l + const x 
|t| 2 / 3 + const' x |t| 4 / 3 ). Casclle et al. [5] also fit the magnetization with a |r| 2 / 3 term. 

Numerical results. We eventually deduce the behaviour of the magnetization 



M.{-\T\)=B\r\ l l 12 {-\n\r\)- 1 



-1/8 



! , 31n(-ln|r| 



4 -In 



T 



-1/8 

(l + a|r| 2 / 3 + ...).(12) 



-l 

with the SS correction or 



Note that the whole bracket corresponds to the correction function of Eq. (3). It is un- 
safe (for numerical purposes) to expand it, since the correction term is not small enough 
in the accessible temperature range |t| ~ 0.05 — 0.25. We have thus to extract an ef- 
fective function J- e ff{— hi l T l) which mimics the real one F{— ln|r|) in the convenient 
temperature range. Defining various effective magnetization amplitudes at different lev- 
els of accuracy, namely B^ f (-\T\) = M_ |r|~ 1/12 (- In It]) 1 / 8 with the CNS leading log- 
term, B%{-\r\) = M_|r|-Vi2 ( _i nM) i/8 (i _ l!±JM) 

B^X(— \t\) = M- 1 -r | 1/12 1- — m | T | £(— In Irl)] 1 / 8 with our universal corrections, we arc un- 
able to recover a sensible 

B(l + a\ T \ 2 / 3 + fe|r| 4 / 3 ) (13) 

behaviour. Of course, it is possible to fit the data to any of these expressions in a given range 
of temperatures, but the coefficients a and b thus obtained strongly depend on the temper- 
ature window and this is not acceptable. Improvement is achieved through the following 
type of fit (instead of Eq. (13)) 

b<?,<-M) = b (i + « + ^fejP ) " 8 a + «H 2/S ). (») 



:m ''"-"V -ta|r| T (-ln|r|) 

The function T{— m|r|) in Eq. (12) now takes the approximate expression T{— ln|r|) ~ 

f 1 + ^TrTR + ^FTSFTF^) ' What is remarkablc is thc stability of the fit to Eq. (14). 
Analysing MC data, we obtain (fit a) C\ = -0.757(1) and C 2 = -0.522(11) which yields 
an amplitude B = 1.1570(1). It is also possible to try a simpler choice in thc narrow 
temperature window, fixing Ci = and approximating thc whole series by the Ci— term 
only (now C\ = —0.88(5), called fit b), which then leads to a very close magnetization 
amplitude B = 1.1559(12). An analysis of SE data gives very similar results. By the way, 
in the case of the magnetization the coefficient b is found to be almost zero and we did 
not include it in Eq. (14) [21]. Note that these estimates follow from a coherent analysis 
of both MC data and SE extrapolations. The errors reported are the standard deviations 
resulting from the fits, since our definition of the temperature window is such that there is 
no finite-size-effect in the r— range considered. For MC data we perform weighted fits (i.e. 
each point is weighted with the inverse statistical error of thc point) while thc fits of thc 
SE data are unweighted. Our major improvement (compared to previous references using 
MC and/or SE data) is not in the quality of thc data themselves, but in thc functional form 
of fit employed which incorporates in an effective function (dependent on the temperature 
window) the effect of the non-universal part of the series of log terms, all universal terms 
being explicitly taken into account. Nevertheless, another source of error comes from the 
effective function itself. We can estimate this additional error in thc following way: we 
compute the amplitudes when changing the coefficients C\ and C*2 of an amount as large as 



x It was proposed in Ref. [16] that ipo = in the Baxter- Wu model and there are no log-corrections. Later 
Kinzel et al. [18] gave supporting considerations. 
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half to twice their optimal values reported above, leading to two estimates for the amplitudes. 
We arbitrarily define the difference of these estimates as the additional error. T_ is the 
only amplitude for which this error is significant. 

We thus obtain a closed expression for the dominant logarithmic corrections which is 
more suitable than previously proposed forms to describe the temperature range accessible 
in a numerical study: 

Obs.(±|T|) ~ Ampl. x x [£(-ln|"r|).F(-hi|<r|)] # 

x (1 + Corr. terms) + Backgr. terms, (15) 

Corr. terms = a|r| 2/3 + 6|t| 4/3 + . . . , (16) 

Backgr. terms = D + £>i|r| + . . . (17) 



where < and # are exponents which depend on the observable considered, and take the 
values 1/12 and —1/8, respectively, in the case of the magnetization. The dots represent 
higher order terms which theoretically do exist, but practically do not need to be included. 

The susceptibility (see Fig. 1) and the energy density can also be fitted to the ex- 
pression above. Our results are summarized in table 1. The efficiency of the fits re- 
lies on the asymptotic form Eq. (15) which in our opinion is based on sufficiently safe 
theoretical grounds. Its validity can furthermore easily be checked (indirectly) through 
the computation of effective amplitude ratios for which all logarithmic corrections have 
to cancel. A specific example is given by the leading behaviour of the energy density 
ratio. The values E(f3) and E(f3*) of the internal energy at dual temperatures are re- 
lated through (1 - e- fj ) E((3) + (l - e~ r ) E((3*) = -2. Defining the quantity A+/A- = 
{E{f3) - E )T a ~ 1 /(E - £:(/3*))(r*) Q ~ 1 , the constant E Q being the value of the energy at the 
transition temperature [22], Eq = E([3 C ) = — 1 — 1/-^, we may expand close to the transi- 
tion point A+(t)/A-(t*) = 1 + (3 - a)a q T + 0(r 1+a ) with a q = -E^ c e~^ = ' n(1 ^ ) . 
This relation, checked numerically, shows that the leading corrections to scaling vanish. 

The universal combinations of amplitudes follow from the results listed in table 1 and 
are summarized in table 2. Fits a and b in these tables refer to the two possible choices for 
the constants C\ and Ci in T{— In |r|) as explained above. 

Conclusion. The main outcome of this work is the surprisingly high values of the ratios 
r + /r_, Tt/T- and clearly far above the predictions of Delfino and Cardy. Note that 
our results are also supported by a direct extrapolation of effective amplitude ratios for which 
most of the corrections to scaling disappear. In the case of the conflicting quantities, this 
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Table 1: Critical amplitudes in the 4-state Potts model. The amplitudes reported correspond to 
the estimates which follow from the analysis of MC data and of SE data with both types of fits. 
The second figures in parenthesis for F_ refer to the additional error discussed in the text. 



F fit # 


B 


r+ 


r_ 


r T 


MC a 


1.1570(1) 


0.03144(15) 


0.00454(2)(20) 


0.00076(1) 


MC b 


1.1559(12) 


0.03178(30) 


0.00484(3)(5) 


0.00073(1) 


SE a 


1.1575(1) 


0.03041(1) 


0.00483(1)(20) 


0.00073(1) 


SE b 


1.1575(1) 


0.03039(1) 


0.00493(1)(5) 


0.00073(1) 



technique leads to r+/r_ = 6.6(3) and 6.5(1), and r T /r_ = 0.160(8) and 0.152(2), using 
respectively fits a and b to fit the MC data. The corresponding figures resulting from fits of 
SE data arc r+/r_ = 6.30(1) and 6.16(1), and r T /r_ = 0.151(3) and 0.148(3). Note that 
the additional source of error is not taken into account in these estimates. 

We believe that our fitting procedure is reliable, and since the disagreement with theoret- 
ical calculations can hardly be resolved, we suspect that the discrepancy might be attributed 
to the assumptions made in Ref. [2] in order to predict the susceptibility ratios. Even more 
puzzling is the fact that Dclfino and Cardy argue in favour of a higher robustness of their 
results for Tt/^- than for r + /r_, but the disagreement is indisputable in both cases. 

Finally, in favour of our results, one may mention a work of W. Janke and one of us (LNS) 
on the amplitude ratios in the Baxter- Wu model (in the 4-state Potts model universality 
class), according to which r+/r_ ~ 6.9 [23]. These results, obtained from an analysis of 
MC data show a similar discrepancy with Dclfino and Cardy's results and a further analysis 
still seems to be necessary. 

* * * 
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A+/A. 


r+/r- 
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source 
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